LANCS-TH/9516 (1995), e-Print Archive: |hep-ph/ 9506336 
IS ELECTROWEAK BARYOGENESIS CLASSICAL? * 



TOMISLAV PROKOPEC tt 

Lancaster University, Physics Department 
Lancaster, LAI 4YB, England, UK 



in 

On 
<3\ 



ABSTRACT 



In this lecture first I present a brief genesis of the ideas on the electroweak baryo- 
genesis and then I focus on a mechanism in which the source of CP violation is 
a C .P-violating field condensate which could occur, for example, in multi-higgs 
—} , extensions of the Standard Model. In the limit of a thick bubble wall one finds a 

V/~j ' classical force acting on particles proportional to the mass squared and the CP 

violating phase. One can study this force in the fluid approximation in which the 
effects of transport and particle decays can be taken into account. A novelty in this 
talk is generalization of the problem to the relativistic velocity. There is a regime 
' in which the final formula for the baryon asymmetry has a rather simple form. 
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)0 1. Prelude 

O 
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Q\ ■ In this lecture I will present some of the recent developments in understanding 

baryogenesis at the electroweak phase transition. 
Qh! I will focus on the work I have done in collaboration with Michael Joyce and Neil 

o ' Turok. The main idea is rather simple: 

<D ' 

In the semiclassical approximation baryogenesis can be described by a CP violating 
classical field condensate. The effect of this condensate can be studied in the fluid 
approximation which takes account of both transport and particle decays. 

s; 

I will also present some of the hystorical ideas that 'seeded' our work as well some 
of the recent developments that I find related or interesting. (I cannot hope for a lack 
of personal bias in my choice.) 

2. Sakharov's Conditions and the Electroweak Theory 

The baryon-to-entropy ratio of the Universe one would like to explain is 

(4-7)lCr n (1) 
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This number is a consistency constraint for nucleosynthesis calculations be in agree- 
ment with the observed primordial elements' abundances. For a pedagogical review 
see for recent developments 0. 

'Before Sakharov' baryogenesis was a fine tuning problem of the intial conditions 
of the Universe. Sakharov realized H that baryons can be generated dynamically 
provided the following conditions are satisfied 

•> Condition 1: there exist processes that violate baryon number 

•> Condition 2: system is out of thermal equilibrium 

•> Condition 3: there exist both C and CP violating processes 

The necessity of Condition 1 is obvious. 

An elegant proof of Condition 2 can be found in @. In thermal equilibrium the 
hamiltonian Ti, is CPT invariant, and the baryon number B is odd under CPT 



(CPT) U (CPT)~ L = H , (CPT) B (CPT) 
In thermal bath baryon number is given by a thermal average 
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We can now insert (CPT)(CPT) 1 = 1 into the trace and use Eq. (g) to get 
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so (B) = 0. 

Next I present an argument which indicates the necessity of CP violation, as 
stated by Condition 3. (A similar proof as above applies, but I want to avoid the 
use of thermal equilibrium.) Recall since the net baryon number is defined as the 
difference between the number of baryons and antibaryons B = b — b and under CP: 
b ^b and b — > b, we have 

(CP) B (CP)' 1 = -B (5) 

In case there is no CP violation the processes which create net b occur with the same 
rate as the processes that create net b and no net B will be created. 

Are the Sakharov conditions realized in the Standard Model? 

We start with the baryon number violation. The Weinberg-Salam theory has on 
the quantum level a remarkable axial current anomaly through which baryon number 
current is violated 
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where Np = 3 is the number of families, and B^ u and and B^ denote the 
SU(2) L and U(1) Y field strengths and fields, respectively, and = ^ upCT W^ 

denotes the dual field strenght. The fact that d^K^ and <9 M fc M are total derivatives 
does not mean that this anomalous current violation can be gauged away without 
any physical consequences by merely shifting the current. The resolution is in the 
nontrivial topological structure of the theory. The integral of Eq. defines how are 
the topology changes in the gauge fields related to the changes of the baryon number: 

2 

AB = N F An cs , n cs = - J K^alS^ (7) 

At zero temperature quantum tunnelling can generate baryons i but the rate is tiny 
~ £-471-/0^ _ j-j. wou i(j no t create even one baryon in all of the (visible) Universe! 

The discovery of the sphaleron Q, a classical solution to the dimensionally reduced 
action which carries a half topological charge and as such, it was argued, may mediate 
baryon number violating processes, lead to a belief that the finite temperature baryon 
number violating rate may be unsuppressed which would make the electroweak 
baryogenesis possible. This was in fact proven by Arnold and McLerran in □ where 
they showed that in the unbroken phase the number of sphaleron transitions (per unit 
volume and time) reads 
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where A(X/ g^,) ~ 1.5 — 2.7; in fact for physically interesting range of A/ to a good 
approximation A ~ 1.7. I have not bothered to quote all pre-exponential factors in 
this relationship; they can be found in 1. The main point is the sphaleron rate is 
exponentially suppressed by the Boltzmann factor with the sphaleron free energy in 
the exponent, which is defined as the saddle point of the dimensionally reduced three 
dimensional action. 

An important condition which ought to be satisfied in order for any electroweak 
baryogenesis model be viable is 

•> Sphaleron erasure: The baryons generated at the transition must not be washed 
out by the subsequent sphaleron transitions in the broken phase §. 

Roughly speaking this means that the sphaleron rate per unit time T sp h in the 
broken phase ought to be smaller than the expansion rate of the Universe H. More 
precisely, since at the time of the phase transition completion, the higgs expectation 
value 0o = ((f)) still grows quite rapidly on the expansion time scale, sphalerons are 
active only a small portion of the expansion time after the completion t3. With this 
in mind one gets slightly less stringent bound on the sphaleron energy 

E sph > 35T (9) 



Taking account of M w = g w <fro/2 and the dependence of the higgs expectation value 
on the couplings (which is roughly (f> ~ g 3 T/Xr) one gets an upper bound on the 



higgs mass. It turns out that with the two loop effetive potential this leads to a 
bound: m H < 30GeV '. If one allows for nonperturbative effects this bound becomes 
less stringent. Recall that the current constraint on the higgs mass is m# > 60GeV 
for the Standard Model higgs. Fortunately this does not mean the end of electroweak 
baryogenesis. Since we are interested in an extended higgs sector which contains CP 
violation and the constraints on the mass of the lightest higgs are much less stringent, 
and the above calculation should be altered accordingly, I believe that at the moment 
no serious discrepancy exists. 

An important calculation of the 'sphaleron' rate in the unbroken phase is still 
missing due to the ill understanding of the infrared sector of gauge theories. Indeed it 
is believed that magnetic fields, which are involved in the processes, are screened at 
the scale of the infamous magnetic mass: m mag ~ g^T, which specifies the 'magnetic' 
field correlation length £ ~ 1/g^T. 

We now present a simple estimate of the rate. There will be instantons of all sizes 
p G {0,£} which represent tunneling trajectories in the configuration space with the 
barrier height of order l/a w p. The smallest barrier is for large instantons p ~ £ with 
the energy ~ T so that there is no exponential suppression to tunnelling; the rate is 
then determined by the prefactor which is of order £ 4 , hence T sp h/V ~ (a w T) 4 . 

A number of numerical studies supports this naive argument 

= K sph (a w T) 4 , K sph ~ 0.1 — 1 (10) 

There are however problems with the simulation. Ambjorn et al. use the real time 
classical theory which is plagued with ultraviolet divergences; the authors hope that 
the lattice cut-off takes care of that. Also the Gauss contraint is not naturally 
imposed. I will take Eq. fllHD as a guidance to the true rate. 

Since the symmetric phase 'sphaleron' rate Eq. ( [XQ|) is much faster than the broken 
phase rate Eq. (0) an efficient baryogenesis mechanism should generate baryons in the 
symmetric phase. This is indeed possible if one takes account of particle transport. 

Regarding Condition 2 perturbative calculations £3 valid for a rather light higgs 
particle m# < 70GeV indicate that the phase transition is first order and likely 
proceed via bubble nucleation providing the required departure from thermal equi- 
librium. To resolve the problem fully one requires to solve for nonperturbative effects 
which is possible only by using lattice simulation. Recently c3 it has been shown that 
the dimensionally reduced theory is well suited for studying. equilibrium properties of 
the transition in a lattice simulation. A preliminary result E-3 indicates that the phase 
transition becomes second order for the higgs mass above about 80GeV. 

What about Condition 31 C is maximally violated. For example no right handed 
neutrino has been observed. 

The main problem for the Standard Model baryogenesis is a small CP violation. 
A natural measure is the CP violating phase from the Kobayashi-Maskawa matrix 
5cp ~ 10~ 19 , and since it is usually argued that ub oc 5cp it seems without an 
additional CP violation the Standard Model baryogenesis is out! This however may 
not be true. Recently a couple of attemts have been made to resolve this formidable 



problem. Shaposhnikov and Farrar III (inspired by a modest CP violation of the 
neutral kaon) have realized that under certain conditions the thermal reflection off 
the bubble wall may lead to a net baryonic current due to the difference in the tree 
level light quarks' masses. This model, since in its original formulation takes no 
account of the possible loss of quantum coherence in the reflection off the wall, has 
been recently jeopardized ^3, G-l Another notable attempt to solve the problem of a 
small CP violation is by Nasser and Turok ta. Based on a linear response analysis 
they argue that there might be instability to formation of a Z field condensate on the 
bubble wall caused by the chiral charge deposit from the reflected top quarks. This 
condensate as we will see below may generate a net baryon number. The longitudinal 
Z field domains of the opposite sign form indiscriminately on the wall and no net 
baryon number is formed unless there is a nonzero CP violation. This gives an 
advantage to the formation of a particular sign condensate. Diffusion of the one sign 
domains against the other enhances the bare CP violation. If this mechanism works 
(as it may for a very heavy top quark) it would result in oc VScp- 
How can one bias baryon number production? 

If there is for example a term in the free energy (possibly dynamically generated) 
that couples to the baryon current 

AF = a^ B (11) 
then a classical formula of the near equilibrium statistical physics gives 

5F 

n B — —isph^B , ^B — T — O-0 (12) 

where \xb is the chemical potential for the baryon number. (In this light Condition 
2 of Sakharov asserts that in thermal equilibrium fiB = 0.) Indeed baryon number 
production is biased if on average a ^ 0. This reasoning sets the stage for the next 
section. 



3. The ideas of Cohen, Kaplan and Nelson 

There are two ideas of relevance for my work: 
•> spontaneous baryogenesis and 
•> charge transport 

The idea of the spontaneous baryogenesis in its intial form 111 is rather simple. 
Cohen and Kaplan assumed that at some early stage in the Universe CPT symmetry 
of the hamiltonian was temporarily violated by a term of form: 

j^d,9J^^9(n b -n- b ), M>10 13 GeV, T<M (13) 

The implication (— ►) applies for a homogeneous field 9. The ilion field 9 (the name is 
supposed to bear reminiscence of the axion) acts as a chemical potential for baryon 



number which, when the expansion of the Universe is taken ito account, drives 5^0, 
hence the name spontaneous. 

The idea got a new impetus when Dine et al. c2l realized that a similar term 
arises on one-loop level in theories which contain a singlet field {e.g. supersymmetric 
theories). The effective action then acquires a term 

' d,sr B (14) 



which may lead to a spontaneous baryogensis at a first order electroweak phase transi- 
tion. Note that in this scenario the electroweak transition is used to drive the system 
out of equilibrium. A similar term occurs at one loop in two higgs doublet theories 
with CP violation in the higgs sector 0. But then Cohen, Kaplan and Nelson 
realized that in these theories even at the tree level there is a term 



ypd^JZ, j£ = ^ 75 vd (15) 

which couples the fermionic hypercharge current Jp, which is closely related to the 
baryonic current, with a CP violating field 9, the gauge invariant relative higgs phase. 
9 acts as a chemical potential for the baryon number. On an expanding bubble wall 9 
is predominantly of one sign, as given by the classical equations of motion, and hence 
it drives the sphaleron rate in a definite direction creating baryons. As originally 
designed this mechanism is not manifestly mass suppressed, i.e. baryon number does 
not vanish in the limit when fermion masses go to zero. In this limit the chiral 
symmetry is restored and baryon creation must stop. I will present a resolution iof 
this paradox below. 

The idea of charge transport was designed @ in a model with heavy (majorana) 
neutrinos and explicit CP violation in the coupling constants. In this model CP 
violation takes place on the bubble wall. As a result of a CP violating reflection on 
the wall an axial current is injected and transported into the symmetric phase where 
the sphaleron rate is unsuppressed Eq. fllPl). This mechanism takes advantage of a 
large CP violation on the bubble wall and of a fast baryon number violating rate in 
the symmetric phase so it is potentially very efficient. Indeed in the subsequent work 
Cohen, Kaplan and Nelson cB have applied the idea to the top quark and managed 
to generate baryon asymmetry as large as ~ 10~ 6 9cp- Since the injected reflected 
chiral top current is the result of a coherent quantum reflection, it is essential that 
top quarks do not scatter on the bubble wall. Also a significant reflection occurs only 
when the wall is thinner than m^ 1 . Both of these lead to a rather stringent constraint 
on the wall thickness L w < 2/T which would require a very strong first order phase 
transition and this is, according to perturbative calculations of the effective thermal 
higgs potential Ii3, highly unrealististic. 

We have realized c3 that this idea might have a 'brighter future' if aplied to 
the particles that couple weakly to the wall and scatter infrequently. Our favorite 
candidate is the right-handed r-lepton since it has a long scattering lenght and diffuses 



very efficiently. We have shown that with 9cp ~ 1 the baryon number produced may 
be consistent with Eq. ([!]). 



4. Two higgs doublet model and Z condensate 

In order to study the effects of CP violating condensate on the particles' dynamics 
on the wall I will rewrite the Lagrangian in a convenient form so that the interesting 
CP violation is expressed in the form of a CP violating Z field condensate. 

Consider first the higgs kinetic term. I can pick a gauge in which the background 
classical solution reads 

* i= ^Ue«.)' * 2= 7§uy (i6> 

In what follows I will ignore the electromagnetic fields (U(1)em remains unbroken so 
there is no obvious reason for a condensate to appear), the charged W ± gauge bosons 
and the higgs excitations both neutral and charged. Then the higgs kinetic term can 
be diagonalized 
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gZ, = g 2 W*- gi B,, 2 = 0? + <f% , g 2 = gl + g 2 2 (17) 

where I ignored the kinetic terms for <pi and 02- Note that when either of the higgs 
vevs vanish the gauge invariant relative phase 9 loses its meaning, as it should. Both 
9 and Z^ 1 are gauge invariant, and we take them as our definition of the CP violating 
condensate fields. 

Consider now the fermionic part of the Lagrangian. Assume for simplicity that 
only $i couples to the fermions via the Yukawa terms. The phase 9\ can be removed 
by performing an anomaly free rotation proportional to [T 3 — Y + |(5 — L)] on the 
fermions (and also on the Higgs field) at the cost of introducing a new term in the 
fermionic kinetic term 

- ig A Z^ 5 ) * - (18) 

where g& = +g/4 (g 2 — g\ + g 2 ) for up-type quarks and (left-handed) neutrinos, and 
qa = for down-type quarks and charged leptons, and 

g A Z, = g A Zf - ^d,9 (19) 

The advantage of this form is that I have reduced the effect the higgs phase to a pure 
gauge Z field condensate in the presence of a (real) mass m = ycj)\/y2. Since the 



axial symmetry is broken this field produces an interesting physical effect which we 
want to study This formulation is suitable for studying the effect of a CP violating 
condensate in the semiclassical limit. 

An alternative formulation is to study the classical evolution of the CP violating 
phase accross the wall which actually neglects the plasma effects. The two-higgs- 
doublet higgs potential may contain CP violating terms 



A 5 (K\ [${$ 2 ] - v x v 2 cos£ ) + A 6 [It [$i$ 2 ] - V!V 2 sin£ 
so that the classical evolution of the relative higgs phase 
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For a slow moving wall d fl d fl 9 ~ can be neglected and we have for the symmetric 
and broken phases 



tan 9 = — — - tan £ , 
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(22) 



so that only if X e 7^ A5 there is a net change in 9 accross the wall, i.e. CP violation 
and Z condensate. For a small CP violation 
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(A 6 - A 5 )A 6 cos£q 
A| cos 2 £ + A| sin 2 £ 



(23) 



which is an upper limit to change in 9 givenby Eq. fl2"Tp. How this may drive 
baryogenesis was first studied in c3 and then in c3. 



5. A free particle on the bubble wall 

We are now ready to study the motion of a free particle on a bubble wall with a 
Z field condensate, i.e. CP violation. Considering a planar wall stationary so that 
Z" = (Z°,0,0,Z 3 (z)). 

First using Eq. (|18|) we write the Dirac equation (in the Fourier space) 



f- QApls ~ m 



* = 



(24) 



which in the chiral representation (where $ ~ [$?r, can be broken into two 
equations for the two component right-handed and left-handed spinors 



E-g A Z°-a-(P-g A Z) V R + m* L 







E + g A Z° + a-(P + g A Z) 



(25) 



where we used 75^^,1, = ±^h,l- Setting the determinant to zero leads to the disper- 
sion relation which is simple in the wall frame where Z^ 1 = (0, 0, 0, Z(z)) 
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(26) 



S z is the component of the spin in the z direction, measured in the frame in which p± 
vanishes. In the WKB approximation, this dispersion relation accurately describes 
particles as they move across a bubble wall - the local eigenstates in Eq. (p6|) shall form 
the basis of our treatment. The particles we are most interested in for baryogenesis 
are left handed particles (e.g. t^s), and their antiparticles (ic's, which are right- 
handed), because these couple to the chiral anomaly. For large \p z \, these are easily 
identifiable in terms of the eigenstates in Eq. (^). Note that they couple oppositely 
to the Z field. 

The group velocity of a WKB wavepacket is determined from the dispersion rela- 
tion by v z = z = dE/dp z , and energy conservation E = implies p z = —d z E. From 
these it is straightforward to calculate the acceleration 



where E and p± are constants of motion. 

This chiral force provided by the Z field effectively produces a potential well which 
draws an excess chiral charge onto the wall, and leads to a compensating deficit in a 
'diffusion tail' in front of the wall. There is net baryon production because B violation 
is suppressed on the bubble wall. 

6. Fluid equations 

I now seek to describe the particle excitations with dispersion relations Eq. (^) as 
classical fluids. I focus on particles with large \p z \ ~ T >> m for the following reasons: 
they dominate phase space, the WKB approximation is valid, and the dispersion 
relation simplifies so one can identify approximate chiral eigenstates. The S z = +|, 
p z < branch, and the S z = — |, p z > branch constitute one, approximately left 
handed fluid L, and the other two branches an approximately right-handed fluid R. 

The Boltzmann equation is: 



where z = d Pz E and p z = —d z E are calculated from the Hamilton equations and 
C(f) is the collision integral. 

In order to study a system close to equilibrium, fluid approximation is reasonably 
accurate, provided 

•> the scatterings within a fluid are more frequent then between any two fluids, 




(27) 



d t f EE dj + Zdj+p z d p j = -C(f) 




•> the distortions in the momentum space are small so that the momentum ex- 
pansion is well behaved (see below); this is the case e.g. for a thick bubble wall. 
In the plasma frame the fluid Ansatz for the distribution function is 



f= eE+s±1 , 5=-y.+p v 8 v + ..., 5» = (-5T,0,0,v) (29) 

where I have for convenience set T = 1 and truncated the expansion at the vector 
term. Since the dispersion relationship Eq. (|2Bp is obtained in the wall frame we 
transform energy and momentum 

E -»• j w (E + v w p z ) , p z -»• "f w (p z + v w E) , p±^p± (30) 

which gives Eq. (p5| ) in the wall frame 

f = e ^(E+vL)+s ± 1 ' 6 = _/i ~ ^ E + V ^ 6T ~ 7«(P* + v*E)v (31) 

There is a subtlety in deriving the fluid equations. In the zero mass limit there must 
be no observable effect (particle number perturbation) due to a nonzero Z field since 
it is pure gauge. However d t p z = —d z E = —{p z ± g A Z) / ' E)d z (g A Z) would generate a 
chemical potential. This is what Cohen, Kaplan and Nelson indentified as the term 
that drives spontaneous baryogenesis. Dine and Thomas B3 and also Joyce, Turok and 
myself il have pointed out that there must be a mass squared suppression. Therefore 
this term is a pure gauge artefact. The proper (gauge covariant) form of / includes 
the kinetic momentum and a shifted 'physical' chemical potential 

k z =p z ± g A Z , fi = /i ± j w v w g A Z (32) 

where the signs are chosen appropriately. The source term 

d t j w (E + v w k z ) = "i w v w k z = j w v w (d z E ± g A Z) (33) 

vanishes in the zero mass limit as it should. I can now write down the linearized form 
of the l.h.s. of Eq. (EH) as 



(34) 



dtf = f'{jw{E + v w p z )) [i w v w k z - d t n - j w (E + v w p z )d t 5T - ^ w (p z + v w E)d t v 

with d t — > d Pz Ed z ~ ^d z in the wall frame. 

Rather then studying the full fluid equations I will simplify the problem by ne- 
glecting the temperature fluctuation. This is reasonably accurate when the scatter- 
ings that destroy ST are much more efficient than the decay processes that destroy /i. 
In the electroweak plasma this is satisfied to iL good approximation. For a detailed 
treatment of the complete fluid equations see Ea. 

In order to obtain the fluid equations I integrate Eq. (|34] ) in the wall frame over 



/ d 3 p and J d z p^/ w {p z + v w E). I have chosen this Lorentz-like combination of E and 



p z since in this case the collision term has a simple form when evaluated in the 
fluid frame: only the velocity perturbation is damped and both chemical potential 
and temperature perturbations drop out. After a lenghty algebra I arrive at the 
fluid equations for particle minus antiparticle perturbations (// = n(L) — fi(L), v — 
v(L) — v( L )) in the rest frame of the wall 



a 2 3V w j w jl' + Iw-^v' = -r M i \fi\ - r* x [/t + 2j w v w g A Z] 



(35) 



azAlwtf - lwV w v' + F 3 = -T v v 

F 3 = 3a M v w (9AZm 2 )' (l - ^ In i±M (36) 

where [p] denotes the signed sum of chemical potentials for particles participating 
in the reaction, prime denotes d z , a n = (n\( n /2n 2 )[l — l/2 n ~ 1 ] (n > 1), a x = 
(ln(2/m)/27r 2 )[lnm], a 23 = a 2 /a 3 = (C 2 /3Cs)[2/3], a 34 = a 3 /a 4 = (Cs/^) [6/7], 
ai4 = 0,1/0,4 = (In 2/24) [8/7] and ( is the Riemann ^-function; to obtain the con- 
stants Oij for bosons one should simply drop the factors in square brackets. 
There are two sources of baryogenesis in these equations 
•> the classical force term oc (gAZm 2 )' in the second equation and 
•> the hypercharge violating processes r* x [2j w v w gAZ] 

The rates and T v damp the chemical potential and velocity perturbations 
respectively. If one for example considers the top quark baryogenesis, the main tree 
level procesess that contribute to are the helicity flipping top-gluon scattering 
with the emission of a charged higgs particle or a charged W boson and permutations 
of these. The process with the higgs emission (absorbtion) is faster due to its large 
Yukawa coupling. For quarks in addition there is also the anomalous strong sphaleron 
process that flips chirality. The hypercharge violating helicity flipping processes are 
(in/nT) 2 suppressed and therefore rather slow (see In this letter I shall focus on 
studying the effect of the classical force. In order to simplify the discussion I shall 
set all decay rates to zero. This is justified provided during the wall passage a small 
portion of the seeded perturbation decays i.e. 

^«r M i,r; 1 (37) 

"VI 

Also it is important that decays do not cut-off the diffusion tail in front of the wall. 
I will comment more on it below. 

I now set all decay terms to zero Y* x — > 0. Eq. fl35D can be integrated to give 

v = 3a 23 v w fl (38) 

since at +00 both jl = v = 0. Eq. ( |3"6"D can then be written as 

/_ ^3 



Dp! + v w jl = s 



3a 23 T v 



D = D(1-^), D = -^- (39) 



where D is the (plasma frame) diffusion constant. The inclusion of the temperature 
perturbations changes D only slightly; in that case D = 1/3T V . We have calculated 
c3 the contribution of the quark-quark gluon exchange tree-level diagram and found 
D ~ 6/T, a similar contribution is expected from the quark-gluon scattering diagram, 
which would half the diffusion constant e3. 

Note that close to the speed of sound v w ~ v s — 1/ v^3 in this approximation the 
fluid becomes stiff, i.e. transport is suppressed (D/v w << L) and there is only local 
baryogenesis: 

A = — = (40) 

which is suppressed as D/Lv w . This result should be taken with caution since the 
limit D ~ might signify the breakdown of the fluid approximation. In a more 
careful treatment one should include other types of perturbations, most notably 5T. 
If one does so one finds that the perturbation at the source may be transported back 
to distances ~ 1/T V . So generically there is no diffusion in front of the wall, but 
particles do diffuse behind the wall. 

In order to write a general solution to Eq. ( |5UD we write down the solution using 
the Greens function method 



d x G + — $( x ~ y) ? G = 6(x — y) exp 



■jj\x-y) 



fi = jGs' (41) 



for D > 0. This form of the solution is quite useful since using partial integration it 
can be expanded in powers of D/v w L for thick walls, or Lv w /D for thin walls. I leave 
this as an exercise. The Greens function in Eq. (^T|) is the response to a 5-function 
source and as such illustrates transport properties of the full solution Eq. pi]) : in 
this approximation when no particles decay all of the sourced perturbation gets trans- 
ported into the diffusion tail in front. This tail extends to D/v w in front of the source 
which is large in the limit of a small velocity. A particle spends on average particle 

Tdiff ~ ^ (42) 

w 

diffusing in front of the wall. If this time is larger than any od the (zero mass) decay 
rates, then the diffusion tail will be suppressed as T^D / v 2 w and accordingly the source 
for baryogenesis will be reduced. 

In order to solve for the baryon number I may integrate Eq. (|12J) and find that 
it suffices to know the integral of fi in front of the wall (for non-local baryogenesis) 
or on the wall (for local baryogenesis). Integrating Eq. fl5§D and assuming efficient 
transport such that v w L/ D << 1, I find that on the wall , 

s 2 In 2 1 — — In ^ h!L 

£ = 7s = — ^—1wVw9AZm 2 V(v w ) l w V(v w ) = ^ — 1 -y h£ - (43) 

D 3C 3 l-3< 023/034 



In order to obtain the baryon number I as anticipated integrate Eq. ( [121) with all 
appropriate constants to find in the plasma frame 

3 , -r f r s p/j 



, n B = n B = -N C / - m L ) (44) 

v w j w J 2 J V 

where Nq = 3 is the number of collors and T sp h/V is the sphaleron rate per unit 
volume. Note that I can easily obtain by integrating Eqs. ( Pop and ( |3"%| ) that J fi = 
Iwaii A + Jo°° A* = 0- This does not mean that baryon production vanishes because 
the sphaleron rate on the wall Eq. (g) is exponentially suppressed, indeed to a first 
approximation I can neglect the baryon production on the bubble wall to find that 
the net baryon to entropy ratio is 

- = J **-jr- v M - 4^— J dz-^-v(v u ) (45) 

where s = ^-g*T 3 is the entropy density, with « 100 the effective number of 
degrees of freedom, and we have restored the units of temperature. One should not 
forget that this result is applicable only in the limit of efficient transport, i.e. not too 
close to the speed of sound, when D ~ and only local baryogenesis is possible. 

This result is remarkably simple - all dependence on the wall thickness and the 
diffusion constant drops out. The dependence on the wall velocity is rather weak; it is 
given by V(v w ) of Eq. ([43]) . For modest velocities V(v w ) ~ (l — v w ) 3 ^ 2 /( 1— 3^023/034) 
so that baryon production becomes more efficient as the velocity increases; the growth 
stops when the condition for efficients transport gets violated: D ~ v w L. The asym- 
metry is quite large for the top quarks: (m t /T) 2 ~ 1/4, so n B /s ~ 10~ 8 n sph 9 C p where 
9c p characterizes the strength of CP violation. The m 2 dependence in Eq. (fHf) means 
that, at least with standard model-like Yukawa couplings, the top quark dominates 
the effect. 

The calculation of the classical force effect above uses the opposite (WKB) ap- 
proximation to those employed in quantum mechanical reflection calculations (thin 
walls @, @). The classical force calculation is in some respects 'cleaner', because 
the production of chiral charge and its diffusion are treated together. The classical 
force affects particles from all parts of the spectrum, mostly with typical energies 
E ~ T, and with no preferential direction, while the quantum mechanical effect 
comes mainly from particles with a very definite ingoing momentum perpendicular 
to the wall: p z ~ m B (Higgs mass). The quantum result falls off exponentially with 
L as the WKB approximation becomes good. The quantum result also has a v ~ x 
dependence coming from the diffusion time in the medium, which the classical result 
loses because the force term is proportional to v w . 
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